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\0 , Abstract 

m 

CN ' Exploiting the idea that the fast partons of an energetic projec- 

'sj- ■ tile can be treated as sources of colour radiation interpreted as wee 

, partons, it is shown that the recently observed property of extended 

^ I limiting fragmentation implies a scaling law for the rapidity distribu- 

tion of fast partons. This leads to a picture of a self-similar process 
where, for fixed total rapidity Y, the sources merge with probability 
rS I varying as 1/y. 

1. It was shown recently [1] that single particle (pseudo) rapidity distri- 
butions in p — p, d — Au, Cu — Cu and Au — Au collisions at 200 GeV cm. 
energy [2] can be described by various superpositions of contributions from 
the "wounded" constituents of the nucleon: a quark and a diquark. The 
form of the contribution from one wounded constituent W{ri) ^ W{y) was 
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determined. Its characteristic feature is that, although peaked in the direc- 
tion of the constituent, it is not restricted to one hemisphere but extends 
over almost full available phase-space [3]. 

Apart from details inessential for our study, this observation shows that 
particle production can be discussed independently for the projectile and for 
the target. This certainly radically simplifies the problem which, originally, 
may seem to be hopelessly complicated. Moreover, it suggests that the main 
effect of the target on the emission of particles from the projectile is a passive 
one: apparently the role of the target is to define the phase-space available 
for the radiation from the projectilqj. 

The analysis of was performed at one single energy. It is interest- 
ing to investigate if (and how) the resulting picture can accommodate the 
recently discovered property of extended limiting fragmentation [H [5] of the 
(pseudo)rapidity spectra of particles produced in high-energy collisions [5],|6]. 
This asymptotic property says that when plotted vs. C, = Y — y where y is 
the rapidity of produced particle and Y is the rapidity of the projectile, data 
at increasing energies approach a common curve (a straight line to a good 
approximation), until a certain value of which increases with increasing Y. 

In the present paper we discuss this problem in some detail and formulate 
a model of particle emission from a strongly interacting composite projectile 
which satisfies extended limiting fragmentation of the spectra and is able to 
lead to a satisfactory description of data. The model exploits the idea of 
distinguishing between the fast partons treated as the sources of radiation 
and the wee partons considered as a colour field emitted from the sources 
O [8] • The emission is described by the bremsstrahlung mechanism where the 
infrared cut-off, necessary to give a physical meaning to the bremsstrahlung 
radiation, is treated as a random variable. 

By demanding that the observed spectra satisfy extended limiting frag- 
mentation, one can show that the distribution of sources must obey a scaling 
law, which expresses the explicit rule how the distribution of sources varies 
as a function of their rapidity y for a given value of Y. Indeed, the smaller is y 
the less is the number of sources. As a consequence of scaling, the same mech- 
anism describes how the distribution evolves at fixed y for varying energy of 
the collision. The scaling law also suggests that the mechanism responsible 
for this structure is a self-similar process in which the sources merge when 

^Of course the roles of projectile and target are completely symmetric and can be 
exchanged one for the other. 
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their rapidity decreases. 

In the next section the mechanism we propose is explained. The scahng 
law following from extended limiting fragmentation is discussed in Section 
3. The self-similar merging picture is described in Section 4. Our results are 
summarized and commented in the last section. 

2. Starting with the idea formulated in we assume that the wounded 
constituent represents a collection of partonic sources, each one emitting soft 
partons which may be represented by a colour field. Denoting the distribution 
of sources by R{y+; Y), where y+ is the rapidity of the sourc^ and Y is the 
rapidity of the projectile, we obtain for the parton distribution 

W{y; Y) = r dy+ H{y- y+) Y) (1) 



where H{y; y^) is the distribution of a// partons coming from a source located 
at ?/+. The final distribution dN/dy is obtained from the sum of contributions 
of the wounded constituents from both particle and target sides. 

The details of H{y; y^) depend on the mechanism of radiation. Here, for 
simplicity, we take the bremsstrahlung or Weizsaecker- Williams distribution 
in its simplest form: 

H{y;y+) = r dy_ P{y^;y+) e(y-y_)0(y+-y) (2) 
Jo 

i.e. a flat distribution in rapiditj|§, extending from y+ till some (in prin- 
ciple arbitrary) infra-red cut-off value P(?/_;|/_|_) is the distribution of 
probability that, for a source located at y+, the radiation is cut at ?/_. 

Physically, the very existence of the cut-off is the consequence of the 
interaction with the target. Since, however, the target is a rather complicated 
object, one cannot hope to evaluate explicitely the actual distribution of ?/_. 
As a statistical approximation we shall assume that it can take any value 
within the allowed kinematic limits with equal probabilitjfl. This idea leads 
to the explicit formula for the distribution P{y-] y+): 

Piy-;y+) dy^ = — Qiy^)Q{y+-y-) dy. (3) 

y+ 



•^Throughout this paper we use the target rest frame. 

■^The general formula for the distribution of partons following from the bremsstrahlung 
mechanism at high energy is a(l — a:)" + ax{l — with x = e^~^+ [H]- The first term 

represents the radiated partons. The second term describes the radiating (leading) one. 
For a = 1 we obtain the flat distribution 

''We shall discuss some generalization at the end of the paper. 
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where the factor in front provides the proper normahzation of the probabihty: 

r P{y^-y+)dy^ = l. (4) 
Jo 

All in all, we finally obtain 

W{y-Y)=y r R{y^-Y)^ . (5) 

Jy y+ 

This formula allows one to determine the distribution of produced partons 
from that of the sources contained in the wounded constituent. 

3. Let us now introduce a convenient formulation of the constraints im- 
posed by the property of extended limiting fragmentation. Starting with the 
distribution W{y, Y) and its derivatives ^ ^^fc'^"* the target fragmentation 
region ?/ ~ "K, it requires on quite general grounds 

dWjy.Y) ^ d'^Wiy.Y) ^ n f t ^ 9 (R\ 
^ \y=Y - est. ; — 0, for > 2 (6) 

when Y becomes large. Indeed, this reflects the fact that at increasing ener- 
gies, the rapidity distribution approaches a common curve which is a straight 
line. We note that from analyticity the higher derivatives cannot be exactly 
zercEl, but their decrease with Y as expressed by the constraint implies 
that limiting fragmentation region extends further in rapidity Y — y with 
increasing total rapidity. 

A more precise mathematical formulation of these constraints which re- 
produces the experimental features described in the introduction can be sum- 
marized in the formula 

W{y,Y) = {Y-yMs)=Y^{s) (7) 

with $(s = 1) 7^ 0. One easily sees that W{y,Y) given by ([7]) satisfies IQ. 
Moreover, one sees that W{y; Y) scales towards a straight line limit at fixed 
rapidity y when Y —>■ oo. More precisely, the slope of the distribution W{y; Y) 
in the limiting fragmentation region is approximately constant (given by 

^The extended limiting fragmentation is an asymptotic property valid at large Y. In 
particular it does not mean that the spectrum depends only on ^ = Y—y even in a restricted 
region in ^. In other terms, one has to decribe not only the limiting fragmentation, but 
also its breaking at some rapidity increasing with Y. 
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$(s = 1)), with significant deviations showing up at some finite Sq and thus 
at a value of rapidity Y — y sq ■ Y increasing with total rapidity. 

Using our formula ([5]), relating the particle distribution W{y,Y) to the 
fast parton distribution R{y]Y), we shall now show that ([7j) implies that 
R{y; Y) obeys the scahng law 

Riy+;Y) = RiyjY) ^ R{s^) . (8) 

Indeed, (E]) implies 

fl(,.y,._,32%I)M = (9) 

oy ds 

so that (IHl) is satisfied. 

The scaling law ([H]) determines the energy dependence of the distribution 
of sources. It shows that the separation between the sources and the field 
scales with the total rapidity. It allows to express the distribution at the total 
rapidity ujY by the distribution at the total rapidity Y, i.e. the evolution of 
the distribution of sources: 

R{y+- ujY) =R(y±-Y)^R (s+ = . (10) 



uo J \ uoY 

This shows explicitely how the idea of [i7j is realized in our mechanism. 

The extended limiting fragmentation realized by this mechanism is shown 
in Fig.l, where the particle rapidity spectra at various values of Y are plotted 
versus Y — y. One sees that extended limiting fragmentation is indeed well 
satisfied. 

4. The scaling law (jS]) suggests that the source distribution is formed in 
a self-similar cascade. Since R is expected to decrease with decreasing s+, 
this process must correspond to merging of fast sources into slower ones. 

Indeed, introducing new variables t+ = \ogy+, T = logF, we see that ([8]) 
implies the relation 

Riy+;Y) = RiT-U), (11) 

and thus an increment of T is equivalent to a negative shift in t+. This is 
typical of the cascade with equal steps and a constant merging probability 
in t+ = logy+. 

In order to make the self-similar structure more explicit, let us give a 
practical example which in fact will also be physically meaningful. 
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Figure 1: Limited fragmentation of the distribution of produced partons. 
dN/dy is represented as a function of Y—y for three incident energies. Left: 
No recombination, R{s+) oc s']_, with 6 = 6. Right: Recombination included, 
R{s+) oc s^/(l + as^) with a = 3,6 = 10 . Note that here, by definition, 
dN/dy = W{y] Y) + W{Y-y, Y), i.e. the sum of contributions of constituent 
quarks from the projectile and the target. 



We consider a self-similar system of fast partons evolving as a function 
of the variable t with constant merging probability. For completion, we also 
consider the possibility that the decrease at high density has to be somewhat 
tamed if the sources screen each other. Hence, it is possible to write a simple 
equation for their distribution, namely 

^ = b[Rit^)-aR\t^)] , (12) 

where b is the merging probability and a the recombination factor. It is easy 
to write the scaling solution of (IT^ as 



Rpe'^'--^^ _ Ro{y+/Yf 
1 + aRoeb(^+-T) 1 + aRo {y+/Yf 



Riy.^ y) = . , = . L.. • (13) 



This solution has all desired features. Indeed the source distribution goes to 
at rapidity ?/ = 0, decreases by merging when y decreases at fixed Y while, 
at a fixed j/, it decreases with increasing Y . In fact this is the solution which 
was used in Fig {1] and it is compatible with the phenomenological information 
from data, as a first approximation. 
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Finally, let us note that such a cascade is equivalent to a cascade with 
equal steps in rapidity but with the merging probability varying as 



h 



(14) 



y 



It is interesting to observe that this result can be obtained if one takes the 
merging probability in the form 



where do is the average (longitudinal) distance between partons in the pro- 
jectile rest frame and d{ii) is their average distance at the rapidity y. Indeed, 
since the longitudinal density is basically determined by the Lorentz fac- 
tor implied by the boost, one may expect that d{y) = e'^do, leading to 
«(?/) oc 1/y. 

One sees that fllSp resembles the well-known perturbative formula for 
the dependence of the effective colour charge on distance. It should be em- 
phasized, however, that there is also an important difference between them: 
while d{y) and d^ refer to longitudinal distances, in the corresponding pertur- 
bative expression the four- dimensional distances are involved. The relation 
between these two situations remains an interesting problem which, however, 
goes beyond the scope of this paper. 

5. Our results can be summarized as follows. 

(i) We investigated the possibility that the particle production can be 
described by two components: one attached to the projectile and another 
one to the target both extending through most of the available phase-space 
in rapidity. This idea is suggested by the success of the phenomenological 
analysis based on the wounded constituent and wounded nucleon models 
[TOj which showed that such an approach can describe the rapidity spectra 
of various processes involving hadrons and nuclei [H |bJ§. 

(ii) Considering one of these components (attached, say, to the projectile) 
we considered the mechanism of particle production suggested in [7j and [8] . 
The partons present already in the rest frame of the projectile are treated as 
sources of the bremsstrahlung radiation producing the colour field. The dis- 
tribution of sources is determined by two effects: (a) The boost from the rest 

^We note that a separation of this sort at high energies has been found in the study of 
nucleus-nucleus interactions in Pomcron field theory [llj . 



a{y) oc 



1 



(15) 



log [do/d{y)] 
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to the laboratory frame implying a substantial increase of the parton density 
in longitudinal direction, and (b) the process of recombination accompanied 
by a decrease of momentum and of the density. 

(iii) The meaningful physical interpretation of bremsstrahlung radiation 
asks for the introduction of an infrared cut-off. To implement the idea that 
the details of the target should not play an important role, we assume that 
the cut-off is randomly distributed through the rapidity space available for a 
given source. We hope that this may be a way to take into account the com- 
plexity and multiplicity of elementary contributions to the soft interaction 
with the target by using statistical arguments. 

(iv) It turns out that the property of extended limiting fragmentation, 
observed in all high-energy processes involving strongly interacting particles, 
implies a scaling for the rapidity distribution of the sources. The scaling 
law states that this distribution depends only on the ratio y/Y and not on 
y and Y separately. This means that the separation between the sources 
and the field is not fully defined by rapidity (as is the case in [7j) but their 
proportion varies with the ratio y/Y . 

(v) The scaling law ([H]) suggests that the recombination process, respon- 
sible for the distribution of sources, is a self similar cascade. This self-similar 
structure differs from the cascade considered in, e.g., [I3] in few important 
aspects. First, it is a merging process instead of a splitting one. Second, it 
is self-similar in logy rather than in y. Third, the transverse momentum of 
partons (which in [13] corresponds to the inverse dipole size) does not need 
to be large. 

Following comments are in order. 

(a) As seen from ([7]), extended limiting fragmentation is an asymptotic 
{Y oo) property. The rate at which it is realized depends on the shape of 
the function 

/ R{s^)^ = [^{s) - ^{1)] = m . (16) 

i — S Js i — S 

One sees that for R{s-^) varying slowly [8j near s+ = 1 (corresponding to 
recombination with sizeable value of a in (fT2|) for the fast partons) one obtains 
a rather fast approach to the limiting line. This is illustrated in Fig.l. Since 

^The existence of a scaling law related to limited fragmentation has been already sug- 
gested in Ref.[12]. 
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the data indicate that hmiting fragmentation is reached already at relatively 
low energies, one may infer that indeed recombination is at work. 

(b) It is worth noticing that the specific structure of the distribution of 
infra-red cut-offs as given by ([3]) is not necessary to obtain the scaling law 
(IH]). Indeed, one can check that any distribution of cut-offs satisfying the 
scaling condition 

P{y.- y+) dy. = Q {s = y./y+) d~s (17) 

also leads to ([8]). 

(c) Our analysis concentrated on the longitudinal properties of the spec- 
tra. It would be very interesting, of course, to elucidate also the role of 
transverse momental§. The key problem should be a better understanding of 
the process of recombination of sources which may or may not affect their 
transverse momentum distribution. Unfortunately, our present approach, 
being purely phenomenological, does not allow to undertake this task. 

d) It should be emphasized that, despite the similarity of certain aspects 
of our approach with perturbative physics, the phenomena we are describing 
are basically non-perturbative. In particular, the scaling property ffTTj) may 
only imply a small longitudinal distance d{y), see f|T5|) . but not necessary a 
small transverse distance [i.e. large transverse momentum) as for saturation- 
like models ([H]). 

A final comment concerns the theoretical implications of the mechanism 
we propos^. Indeed, this mechanism may give some interesting guidelines 
for an exploration of non perturbative properties. One aspect is a theoretical 
hint for the "superposition" property of target and projectile components 
(see e.g. pQ]). Another one is the observation that a self-similar interaction 
is at the origin the longitudinal structure of hadrons boosted to very high 
momenta, which may represent an important step in understanding soft QCD 
at high energies. 
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®In the discussion of extended limiting fragmentation in the framework of CGC ap- 
proach, transverse momenta play an important role |,14j . 

^See, e.g. Refs.^lSj for some recent theoretical discussions on extended limiting frag- 
mentation. 



9 



References 



A. Bialas and A. Bzdak, Phys. Rev. C 77, 034908 (2008). For a review, 
see A.Bialas, J. Phys. G: Nucl. Part. Phys. 35 044053 (2008). 

PHOBOS Collaboration: B. B. Back et al. , Phys. Rev. C 65, 061901 

(2002) ; ibid. 70, 021902 (2004); ibid. 72, 031901 (2005); ibid. 74, 021901 
(2006); R. Nouicer et al, J. Phys. G 30, S1133 (2004). 

A. Bialas and W. Czyz, Acta Phys. Polon. B 36, 905 (2005). See also: 
M. Gazdzicki and M. Gorenstein, Phys. Lett. B640, 155 (2006). 

J. Benecke, T. T. Chou, C. N. Yang and E. Yen, Phys. Rev. 188 (1969) 
2159. 

PHOBOS Collaboration: B. B. Back et al. , Phys. Rev. Lett. 91, 052303 

(2003) ; Nucl. Phys. A 757 (2005) 28; B.Back et al., Phys. Rev. 74, 
021901 (2006); B.Alver et al., arXiv: 0709.4008. 

UA5 Collaboration: G. J. Alner et al. , Z. Phys. C 33, 1 (1986). 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994). 

A. Bialas and M. Jezabek, Phys. Lett. B 590, 233 (2004). 

L. Stodolsky, Phys. Rev. Lett. 28, 60 (1972). 

A. Bialas, M. Bleszynski and W. Czyz, Nucl. Phys. B 111, 461 (1976); 
A. Bialas, W. Czyz and W. Furmanski, Acta Phys. Polon. B 8, 585 
(1977). 

S. Bondarenko and L. Motyka, Phys. Rev. D 75, 114015 (2007); L. Mo- 
tyka. Acta Phys. Polon. B 37, 3533 (2006). 

L. McLerran, Acta Phys. Polon. B 37, 3237 (2006). 

A. H. Mueller, Nucl. Phys. B 415, 373 (1994); 

E. lancu and R. Venugopalan, "The color glass condensate and high 
energy scattering in QCD", arXiv:hep-ph/0303204; 
R. Peschanski, Acta Phys. Polon. B 37, 3511 (2006). 



10 



[14] F. Gelis, A. M. Stasto and R. Venugopalan, Eur. Phys. J. C 48, 489 
(2006); F. Gelis and R. Venugopalan, Acta Phys. Polon. B 37, 3253 
(2006). 

[15] J. Bias de Deus and J. G. Milhano, Nucl. Phys. A 795, 98 (2007), Phys. 
Lett. B 662, 129 (2008). 

J. Cleymans, J. Strumpfer and L. Turko, "Longitudinal Scaling and the 
Thermal Model," arXiv:0712.2463 [hep-ph]. 



11 



